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Abstract Low and oscillatory wall shear stress is widely assumed to play a key role in the initiation and development of 

atherosclerosis. Indeed, some studies have relied on the low shear theory when developing diagnostic and treatment 
strategies for cardiovascular disease. We wished to ascertain if this consensus is justified by published data. We per- 
formed a systematic review of papers that compare the localization of atherosclerotic lesions with the distribution of 
haemodynamic indicators calculated using computational fluid dynamics. The review showed that although many arti- 
cles claim their results conform to the theory, it has been interpreted in different ways: a range of metrics has been 
used to characterize the distribution of disease, and they have been compared with a range of haemodynamic factors. 
Several studies, including all of those making systematic point-by-point comparisons of shear and disease, failed to find 
the expected relation. The various pre- and post-processing techniques used by different groups have reduced the 
range of shears over which correlations were sought, and in some cases are mutually incompatible. Finally, only a 
subset of the known patterns of disease has been investigated. The evidence for the low/oscillatory shear theory 
Is less robust than commonly assumed. Longitudinal studies starting from the healthy state, or the collection of 
average flow metrics derived from large numbers of healthy vessels, both in conjunction with point-by-point 
comparisons using appropriate statistical techniques, will be necessary to improve our understanding of the relation 
between blood flow and atherogenesis. 
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1. Introduction 

The non-uniform distribution of atherosclerotic lesions in the arterial 
system had been attributed to variation in blood flow characteristics 
as early as the nineteenth century.^ Today, the concept that local 
haemodynamic factors account for the patchiness of the disease has 
become almost universally accepted by researchers in the field.'^ 
Acceptance emerged after several decades during which the nature 
of the flow features that drive atherogenesis was the subject of con- 
siderable debate. While Fry^ suggested that high wall shear stress 
(WSS) leads to endothelial damage and a subsequent pathological 
response, Caro et ai"* proposed that high WSS was protective, 
lesions instead occurring in regions of low WSS. In the 1980s, flow 
reversal was identified as another key factor, and an oscillatory 
shear index (OSI) was introduced by Ku et al.^ to quantify it. This 



index has the time-averaged magnitude of WSS in the denominator, 
so high OSI tends to be associated with low mean shear; the low 
and oscillatory shear theories have become fused to some extent. 

In their 1992 review of fluid dynamics at arterial bifurcations, Lou and 
Yang^ noted, 'The theory of low shear stress prevails at the moment, al- 
though the correlation between high-atherosclerosis probability areas 
and expected regions of low wall shear stress is far from perfect. It is, 
therefore, not prudent to accept one theory over another'. Twenty 
years later, the low and/or oscillatory WSS hypothesis has become the 
consensus mechanism for the initiation of atherosclerosis.^'® Some 
studies even rely on the theory for the development or assessment of 
potential diagnostic or treatment strategies.''^" We wished to ascertain 
if the literature really does justify the consensus. 

Blood flow has been intensively studied in recent years. The tech- 
niques of computational fluid dynamics (CFD) have commonly been 
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employed.'^'""^'' In this approach, numerical methods are used to ap- 
proximate the flow fields in arteries, which are too complex for ana- 
lytic solution; haemodynamic factors such as WSS are derived in turn 
from the flow fields. The increasing powencost ratio of computers 
and the ready availability of techniques for incorporating image- 
derived geometries into CFD models have made this approach 
more reliable than methods based solely on in vivo measurement. 

The aims of the review presented here were to identify articles that 
compared CFD-derived metrics with lesion distributions, to collate 
their conclusions, and to gauge the strength of the evidence for the 
consensus view of the relationship between WSS and atherosclerosis. 
To increase the rigour of the review, the literature search was 
conducted in a systematic fashion. 

2. Methods 

A systematic search was conducted on the PubMed database^^ for articles 
published online before 21 March 2012 whose abstract and/or title con- 
tained at least one l<eyword from each of the following three groups: 

• athero, atheroma, atheromata, atheromatous, atherosclerosis, athero- 
sclerotic, atherogenesis, atherogenic, intima-media thicl<ness, intimal 
thickness, wall thickness; 

• shear; 

• CFD, computational, computation, compute, computed, computing, 
numerical, calculation, calculate, calculated, calculating, simulation, 
simulate, simulated, simulating. 

The articles were then evaluated against a set of inclusion and exclusion cri- 
teria. To be included, compliance with the following criteria was required: 

• the work was an original contribution; 

• CFD was used to simulate blood flow in one or more vascular 
geometries; 

• at least one of these vascular geometries was anatomically realistic — i.e. 
it was based on measurements or images, was sufficiently extended and 
included the features necessary to capture the physiologically relevant 
3D flow dynamics; 

• the localization of atherosclerotic lesions in the region of interest was 
based on the distribution of an established marker of the disease, 
presented either as an original contribution in the article, in a compan- 
ion paper from the same group, or in a cited article from a different 
research group. The species used for the disease localization study 
and CFD analysis had to match; 

• the patterns of shear-related parameters in the region of interest were 
shown or described in detail, and compared with the localization of 
disease in the same region; 

• the relation between blood flow and the initiation of atherosclerosis 
was discussed. 

Excluded from the review were: 

• studies that performed CFD on severely stenosed geometries and 
compared the resulting WSS maps with the spatial distribution of 
occlusion in these geometries. With such advanced disease, it cannot 
be assumed that variations in shear are a cause rather than an effect 
of the lesion; 

• studies in which the geometries for the CFD analysis were obtained 
after angioplasty or endarterectomy in the region of interest, since 
these techniques cannot guarantee an accurate reproduction of the 
original 'healthy' vascular geometry; 

• studies in which a stent, cuff, graft, or other foreign object was intro- 
duced in the region of interest, since contact with the foreign material 
or stresses induced by the presence of the object may trigger mechan- 
isms unrelated to normal atherogenesis. Studies in which the geometry 



of the vascular region of interest was itself surgically altered, e.g. by 
anastomosis, were also discarded since effects related to the surgery 
may act as confounding factors; 

• studies in which the disease was mechanically induced, e.g. by denuda- 
tion of the endothelium. Again, these studies may not replicate the 
normal initiation of the disease; 

• studies in which gene expression, white blood cell density, or low-density 
lipoprotein (LDL) permeability were used as surrogates for lesions, since 
these properties are not established markers of the disease. 

Supplementary material online. Figure SI charts the search criteria and 
explains the number of articles eliminated at each stage. 

3. The selected literature 

Twenty-seven articles^'^'^^'^^"'"' (see Supplementary material online, 
Table S1) were selected from an original set of 406. (Supplementary ma- 
terial online. Table SI also lists eight associated papers^^"^® containing 
lesion data that were not presented in the selected articles themselves.) 
The selected articles originated from 15 different institutions and studied 
human, porcine, rabbit, and murine blood vessels; coronary arteries 
were the most frequently studied vessels, but the carotid bifurcation, 
aorta, and internal common carotid were also investigated {Table 1). 
In the animal studies, interventions were required to initiate or acceler- 
ate the disease process: all animals were cholesterol-fed, mice had an 
apolipoprotein E or low-density lipoprotein receptor (LDLR~'~) null 
mutation, and, in one case, the flow conditions were modified by surgi- 
cal induction of aortic valve regurgitation."^^ 

4. Findings 

Section 1 uncritically presents the results for the association between 
blood flow and atherosclerosis given in the reviewed articles. Subse- 
quent sections also consider the methods that were used. Since differ- 
ent research groups applied different levels of data reduction. Section 
2 investigates how data-processing methodologies differed between 
articles. Given the high sensitivity of blood flow to geometric 
changes, Section 3 details at what stage of the disease process the 
flow was characterized. Finally, Section 4 gives an overview of the dif- 
ferent techniques used to assess disease location. In each of these 
three sections, implications for the evidence concerning the role of 
blood flow in early atherosclerosis are discussed. 

4.1 Many but not all articles interpret their 
results as supporting the low/oscillatory 
shear theory 

4.1.1 Studies of steady flow 

Based on investigations of steady-state flow, in which the time-varying 
character of the flowing blood was not simulated, nine papers 



Table I Species and vascular regions examined In the 
reviewed articles 





Human 


Porcine 


Leporine 


Murine 


Coronary arteries 


12 


4 






Carotid bifurcation 


5 








Aorta 


1 




2 


2 


Internal carotid 


1 
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concluded that their results confirm the theory that low WSS pro- 
motes atherogenesis; 

• Krams et ol.^^ and Wentzel et ol.^^''^^ reported an inverse correl- 
ation between wall thickness and steady WSS; 

• Goubergrits et d}° saw a good visual correlation between loca- 
tions of fibrous plaques and regions of low WSS; 

• Chatzizisis et al?^ found that fibroatheromata and intermediate 
lesions develop in vascular segments with lower baseline WSS 
than segments with minimal lesions; 

• Koskinas et d?^ observed a greater increase in maximum intima- 
medial thickness in vascular segments with low baseline WSS com- 
pared to vascular segments with high baseline WSS; 

• Olgac et d?"^ reported a region of low WSS that coincided with the 
location of a plaque; 

• Siogkas et d}^ noted the largest plaque formation in a region 
of low WSS; 

• Chatzizisis et d?^ found that lipid accumulation was enhanced in 
segments with low baseline WSS. 

In contrast, a quantitative comparison by Gijsen et al?^ did not result 
in an overall relationship between WSS and wall thickness. 

4.1.2 Studies of time-varying flow 

The first papers to report on time-varying flow, by Perktold and 
Resch^^ and Perktold et o(.,^^ concluded that the location of athero- 
sclerotic plaques correlated with the extent and location of a recircu- 
lation zone (and consequent low WSS) occurring during part of the 
cardiac cycle. However, subsequent studies focused on time-integrated 
effects. 

4.1.2.1 Time-averaged WSS 

Three papers concluded that low time-averaged wall shear stress 
(TAWSS; WSS magnitude averaged over the cardiac cycle) co-located 
with atherosclerosis or was associated with its rate of progression: 

• Suo et ai^^ concluded that low TAWSS occurred at sites with an 
increased incidence of lesions; 

• Zhu et o(.^^ compared TAWSS magnitudes in two mouse strains 
and found lower values in the strain having greater athero- 
susceptibility; 

• Samady et al.^^ found that plaque area increased over time 
in vascular segments with low TAWSS, but decreased in 
intermediate- and high-TAWSS segments. 

However, two studies did not support the low TAWSS theory: Joshi 
et ol.^^ could not find a significant correlation between intimal thick- 
ness in three out of four arteries, and van der Giessen et ai^^ could 
not establish a correlation with plaque thickness in a mildly diseased 
artery. 

4.1.2.2 Other sliear-derived metrics 

The OSl, first proposed^ to characterize flow pulsatility and later gen- 
eralized to three-dimensional flow,'" was computed alongside 
TAWSS in many of the reviewed papers. 

The conclusions for TAWSS and OSl were in agreement in four 
papers that found an association with early atherosclerosis: Taylor 
et o(.^^ and Buchanan et ai^^'^^ noted that regions of low TAWSS 
and high OSl both coincide with a high probability-of-occurrence of 
early atherosclerotic lesions, and Olgac et al.^° reported that a 
region of low TAWSS and high OSl coincided with the location of 
a plaque. 



One paper found agreement between TAWSS and OSl but no 
association with atherosclerosis: Steinman et ol.^' could not find a 
quantitative general relationship between TAWSS or OSl and wall 
thickness. 

In other cases, the conclusions for TAWSS, OSl, and other metrics 
disagreed, or at least showed discrepancies: 

• Augst et ol.^* computed TAWSS, OSl, and wall shear stress angle 
gradients (WSSAG). They reported a significant relationship 
between WSS or WSSAG and intima-medial thickness in one 
part of the carotid but not another, and no relation with OSl 
in either region; 

• Hoi et al.^^ concluded that OSl and the relative residence time 
(RRT) could explain changes in plaque distribution as a result of 
changes in flow conditions, whereas TAWSS could not; 

• According to Zhang et ol.,^' an OSl-based risk factor was a better 
predictor than a TAWSS-based one of subsequent stenosis; 

• Knight et ol.^^ and Rikhtegar et d."*" found that the largest number 
of plaque locations could be predicted using a low TAWSS thresh- 
old, but that OSl and RRT produced fewer false negatives. 

4.1.3 Summai-y 

Table 2 provides an overview of the conclusions presented in the 
selected articles. Most stated that their results were in agreement 
with the low/oscillatory shear theory, but five^'''^^''^^''^^'^'^ did not 
find a significant overall relationship. 

4.2 Different and sometimes conflicting 
levels of data reduction were required 
to reach a conclusion 

4.2.1 Descriptive analysis 

In the two earliest papers that satisfied all selection criteria,^^'^* the 
comparison of flow and disease was essentially descriptive; differences 
in plaque location in carotids with different bifurcation angles^^ could 
be explained by differences in the extent and location of a recircula- 
tion zone. Such analysis does not provide any statistical evidence for 
the relation between WSS and atherogenesis, in part because of the 
qualitative nature of the clinical observations on which these studies 
were based. (The clinical observations were quantified later^°). 

4.2.2 Visual comparison of maps 

Eleven articles^^'^®''^"''^'^''^^''^®"^'^'^^ presented a visual comparison 
between maps of flow-related parameters and disease indicators. In 
all cases, this approach led to a conclusion in favour of the low/oscil- 
latory shear theory. However, anomalies that run counter to this 
trend are readily visible in the maps presented in at least three of 



Table 2 Overview of the conclusions presented by the 27 
reviewed articles 



Theory 


Supported 


Not supported 


Low WSS (steady flow) 


9 


2 


Low instantaneous WSS 


2 


0 


Low TAWSS 


3 


1 


TAWSS and OSl 


4 


2 


TAWSS or OSl, or RRT 


4 


0 
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the articles.^ ^'■^"'■^■^ Additionally, although Hoi et al?^ advocated the 
use of RRT as a surrogate marker of atherosclerosis, they conceded 
that this metric failed to predict the absence of plaques in the prox- 
imal descending thoracic aorta of LDLR~'~ mice. Another of these 
studies'^^ inadvertently transposed the labelling of branches in one 
map but not in the other; the degree of correlation needs re- 
examination. Thus, the weight of evidence obtained v^ith this 
method is open to question in at least 45% of the studies. 

4.2.3 Point-by-point comparison 

Quantitative point-by-point comparison is a more rigorous form of 
analysis. Steinman et o/.,^' Joshi et ai,^^ Wentzel et o(.,^^ Gijsen 
et al.,^^ and Augst et ai}^ performed point-v\/ise comparisons of 
WSS with intimal, intima-medial, or wall thickness. Spatial resolution 
varied from summary measures for each 1.5-mm wide rectangular 
patch'^^ to one data point per node on the surface of the finite 
element mesh.^' Importantly, none of these studies found a significant 
correlation. 

4.2.4 Axial averaging 

A study by Wentzel et at}'* was one of those in which a 
point-by-point comparison did not yield the inverse relationship 
between wall thickness and WSS expected according to the consen- 
sus view. The authors went on to average the data in the axial direc- 
tion.^^'^^ They argued that this approach should be used in studies 
(such as theirs) where WSS is computed for the geometry of the dis- 
eased vessel but used to explain the initiation of disease, since it will 
reduce the effect of variations in luminal shape or size which would 
not have occurred in the disease-free vessel. An inverse relation 
between WSS and disease was obtained. 

This technique also resulted in a significant inverse WSS-intimal 
thickness relation for Krams et ol.^^ but not for van der Giessen 
et al?'^ In a related technique, Joshi et alP eliminated axial variation 
by normalizing the data per cross-section; as already noted, they 
found a significant inverse WSS-intimal thickness correlation in only 
one out of four cases. 

4.2.5 Circumferential averaging 

Although Wentzel et aiP advised against the use of circumferential 
averaging, since it obscures curvature-induced circumferential varia- 
tions in WSS, Chatzizisis et ai.}^'^^ Koskinas et al.}^ and Samady 
et ai?^ obtained a significant inverse correlation using this technique. 

4.2.6 Selective analysis 

A number of articles limited their analysis to a subset of the data. The 
pre-selection was based either on geometrical features or on values of 
shear- or disease-metrics. 

• Wentzel et al}^'^'* Chatzizisis et o(.,^^'^' and Koskinas et ai?^ 
excluded data from regions of branching; 

• Zhang et a(.^' visually preselected locations of low TAWSS and high 
OSI, and then performed a more quantitative analysis on the 
subset; 

• Chatzizisis et ai}^'^^ and Koskinas et al?^ analysed only segments 
with uniform intima-medial thickness and WSS. 

In these cases, the pre-selection will have restricted variation in shear 
stress or disease severity and may, therefore, have obscured import- 
ant aspects of the relation between WSS and disease. 



4.2.7 Thresholding 

Although most articles did not specify a threshold below or above 
which a haemodynamic factor was believed to promote atheroscler- 
otic disease, five articles presented their haemodynamic results in 
such a categorical sense. 

• Koskinas et al.}^ who studied the flow in porcine coronary arteries, 
fixed the WSS threshold at 1.2 Pa, a value previously used in a 
study of flow in human coronaries;^^ 

• In a second study of porcine coronary arteries by the same 
group,^^ a threshold of 1 Pa was chosen; 

• Samady et ai.}^ in a study of human coronary arteries, defined low 
TAWSS as <1 Pa and high TAWSS as >2.5 Pa. 

The significance of the results obtained in such studies will obviously 
depend on the choice of threshold. Although it is commonly assumed, 
as in these studies, that normal arterial shear lies in the 1.0-1.5 Pa 
range, allometric arguments demonstrate that this is not a universal 
property.^'^ A better approach might be that of Knight et ai?'* and 
Rikhtegar et o/.^° who dichotomized haemodynamic variables by 
balancing the number of correctly predicted plaque locations and 
false positives for 10 patients. 

4.2.8 Summary 

figure 1 summarizes the various data-processing methods employed in 
the reviewed articles. Straightforward point-by-point comparisons did 
not yield results consistent with the consensus low shear theory. 
Although averaging did produce some significant results, averaging 
techniques advised against by some groups were essential for 
others' analyses. Pre-selection of data may have restricted the range 
of shears and disease severities that were considered and hence 
limited the validity of the conclusions. Thresholding of haemodynamic 
parameters needs to be based on objective criteria; furthermore, the 
categorical approach and the continuum approach represent different 
interpretations of the low shear theory — only mechanistic studies can 
determine which (if either) is correct. 

4.3 Blood flow was assessed at different 
stages in the disease process 

Small changes in geometry can induce significant alterations of 
WSS.^^'^^ The validity of studies which compute patterns of shear 
therefore depends on the anatomical fidelity of the reconstructed vas- 
cular geometry. 

4.3.1 Subject-specific studies 

4.3. 1. 1 Single geometry per subject 

Studies of human vessels have often used the same geometry for the 
assessment of disease and for the computation of WSS. Obviously, 
lesions have to be present in the segment in order to allow such an 
investigation. In some cases, the geometry was obtained from indivi- 
duals who were free of cardiovascular symptoms or, post-mortem, 
from individuals whose primary cause of death was not cardiac 
related.^ '■■^■^'■^^ On the other hand, many studies used subjects with 
known atherosclerosis, such as patients who went to clinic with a 
lesion requiring evaluation or treatment.^^''^°''^^''^''''^^'"^'^ 

In studies of this type, there is an explicit or implicit reliance on 
compensatory remodelling^^'^^ to ensure that the presence of 
lesions did not alter the disease-free geometry. In the absence of 
compensatory remodelling, this approach could obscure an inverse 
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Figure I Schematic overview of the various data-processing methods employed in the reviewed articles, illustrated on lesion and WSS maps taken 
from Cremers et a(.'^ and Peiffer et ol.," respectively. The summarizing table on the right indicates that the correlation between the lesion and haemo- 
dynamic distributions, as perceived by the authors of the reviewed articles, increased for increasing levels of data reduction. Note that articles were 
double counted if multiple analysis techniques were reported. 



correlation between early atherosclerosis and WSS since a narrowed 
lumen leads to an increase in velocity (for the same flow rate). 

4.3.1.2 Virtual plaque removal 

The concept of removing plaques virtually from diseased geometries 
is attractive. Olgac et ol.,^''^° Knight et ol.,^"* and Rikhtegar et ai'*° used 
an image-based approach for virtual plaque removal, but acknowl- 
edged the potential for under- or overestimating the size of the 
healthy lumen. 

4.3.1.3 Longitudinal studies 

Samady et al.^^ were the only authors in the final selection to report a 
longitudinal study in human subjects: they re-scanned each person 6 
months after the first examination. However, this study focused on 
progression rather than initiation of disease: the baseline geometries 
were taken from patients rather than healthy volunteers. Longitudinal 
studies are easier in animals, but again the focus has been on progres- 
sion rather than initiation.'^^'^^'"^^ 

4.3.2 Multiple group studies 

A different option is to abandon the subject-specific approach. Instead, 
the distribution of lesions is assessed in a group of diseased vessels, 
blood flow is simulated in geometries obtained from a group of 
healthy vessels, and the comparison is performed between group 
means. Each group needs to be big enough to make inter-subject vari- 
ation unimportant. In most such studies, however, the haemodynamic 
group has been limited to just one geometry;^'^'^'*'^^'^®''^'^''^^'^^ only 
Zhu et ol.,^^ using mice, performed a full statistical comparison. 



4.3.3 Summary 

There have been no longitudinal studies that examined the initiation 
of disease in people. Instead, researchers have focused on progres- 
sion, or have made assumptions about the nature of vascular remod- 
elling, about how to remove lesions virtually, or about inter-subject 
variation. The high sensitivity of WSS stress to geometry makes the 
validity of these assumptions critical. All five studies that failed to 
find a relation between flow metrics and early lesions used the dis- 
eased geometry for flow mapping. However, not all studies using 
this technique failed to find a relation. 

4.4 Disease localization was assessed 
at different stages in the disease process 

Atherosclerosis was mapped through the identification of one or 
more of its constituents or by measuring geometric anomalies in 
wall and/or luminal dimensions. Some studies investigated the likeli- 
hood that lesions will occur at a particular level of WSS, whereas 
others investigated the severity of the disease in relation to shear. 

4.4.1 Disease assessment based on lesion composition 

4.4.1.1 Lipid staining 

Recognizing that the onset of atherosclerosis is characterized by lipid 
accumulation in the inner wall,^^ a number of research groups made 
use of lipid staining with oil red O or Sudan IV to identify the earliest 
stages of the disease. This technique was applied in the studies of 
Taylor et ol.,^'* Goubergrits et al.,^° Suo et aL,^^ Chatzizisis et d.,^^'^^ 
and Hoi et o(.,^^ or their companion studies. 
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4.4.1.2 Calcium scoring 

As the disease progresses, other plaque components appear. Olgac 
et a(.,^''"^° Knight et al.,"^^ and Rikhtegar et aL'*° focused on the calcium 
content of the wall, which they obtained by computed tomography. Dis- 
advantages ofthis method are that it cannot locate non-calcified plaques 
and that the calcified areas may not exactly match the locations where 
the disease was initiated if there is a continuous two-way interaction 
between flow and the development of the disease. In addition, calcifica- 
tion without atherosclerosis is possible; confirmation of plaques by an 
experienced clinician, as was carried out in the studies by Knight 
et d.^^ and Rikhtegar et al.^° is, therefore, important. 

4.4.2 Disease assessment based on geometric changes 

Early wall thickening might not be detectable from an assessment of 
luminal dimensions, due to compensatory remodelling. Hence, a 
large number of the selected articles opted instead to use intimal, 
intima-medial, or wall thickness as a measure of disease. 

Joshi et al.,^^ Zhu et ol.,^^ and Chatzizisis et ol.,^^ or their companion 
papers made use of post-mortem preparations to assess thickness. 
Imaging techniques, such as magnetic resonance imaging (MRI),^''"^' 
ultrasound,'^^ and intravascular ultrasound, ^^'■^^'■^^'■^^'■^^'"^■^'^^■^^■^^ have 
been used for the same purpose in vivo. There are disadvantages, 
however. With intima-medial and wall thickness, it is not possible 
to distinguish between intimal thickening due to atherosclerosis, 
medial thickening in response to elevated blood pressure^^ or adap- 
tive intimal thickening unrelated to atherosclerosis; medial thinning 
under lesions may also cause problems. Imaging methods that 
combine the measurement of intimal thickness with the identification 
of intimal composition in vivo would be ideal. 

4.4.3 Probabllity-of-occurrence vs. plaque severity 

Studies of the relation between flow and disease have variously inves- 
tigated whether or not atherosclerotic lesions occur in certain 
haemodynamic environments,^®''^°''^'^''^''^°'^^'"^^ or how likely it is to 
find lesions in these regions,^'*''^^'"^^'^' or have interpreted the low/os- 
cillatory shear theory in a dose-dependent way by asking whether 
disease is more advanced in these regions.^^'^'"'^^''^^"^'''^^"^^'^*'^'' 

For completeness, we note that an additional category of articles 
have addressed a slightly different question by asking if certain haemo- 
dynamic conditions can protect against the disease.^^'^° 

4.4.4 Summary 

A variety of techniques were used to assess the location of atheroscler- 
otic disease. Whereas lipid staining will have picked up early signs of ath- 
erosclerosis, this is not the case for approaches involving calcium 
scoring or luminal narrowing. Moreover, measurements of intimal, 
intima-medial, and wall thickness cannot specifically identify intimal 
thickening due to atherosclerosis; the five articles which did not find a 
correlation between WSS and atherogenesis used these methods, al- 
though some articles using these methods''^'^'''^'''^'''^^'^^'^*"^'' did find 
a relation. Finally, it is not clear if the low/oscillatory shear theory 
should be interpreted in terms of the frequency or severity of disease, 
or whether both measures are satisfactory. 

5. Discussion 

The low shear stress theory was heterodox when first proposed. It 
challenged the prevailing high shear stress theory, which appeared 
to account for the pattern of disease observed in various non-human 



species and which was associated with the intuitively appealing mech- 
anistic hypothesis that high shear stress damages the endothelium, 
allowing excessive entry of plasma lipoproteins into the wall. 

The involvement of oscillatory shear emerged from the studies of Ku 
etol.,^ who examined lesion distributions in post-mortem human arter- 
ies and used laser Doppler anemometry to measure velocity profiles in 
anatomically accurate models of those arteries. Indeed, the relation 
between intimal thickness and OSI was much stronger than that 
between intimal thickness and mean shear. Today the low/oscillatory 
shear theory is widely regarded as established, and it is contrary views 
that are regarded as heterodox. 

The present systematic review considers whether the evidence jus- 
tifies that consensus. It focused on studies that employed CFD to 
assess patterns of wall shear because these are currently the most re- 
liable methods. (Recent improvements in the accuracy of flow meas- 
urement by MRI^® and other techniques^' may change this in the near 
future.) 

We cannot exclude having overlooked some relevant papers, 
despite the rigorous methodology. Nevertheless, reliable conclusions 
can be drawn from those that were reviewed. First, the majority of the 
studies supported the low/oscillatory shear theory but some did not. 
Secondly, some found the same relation with disease for low WSS 
and OSI, but others did not. Thirdly, close inspection shows that 
the theory has been interpreted and tested in different ways. This 
third point could be regarded as evidence that the theory is robust 
and insensitive to precise formulation — supporting evidence was re- 
peatedly obtained despite the following methodological issues: a 
variety of metrics were used to represent the disease distribution; 
blood flow was assessed at different stages in the disease process; 
data reduction techniques limited the range of shears or disease se- 
verities under consideration, or had mutually contradictory rationales; 
and WSS was treated as a continuous or a categorical variable. An al- 
ternative reading is that there is a lack of unanimity between studies, 
and that the diversity of methods weakens rather than strengthens 
support for a single, precise hypothesis. According to this viewpoint, 
the low/oscillatory theory may very well be true but — taken to- 
gether — the evidence for it is less clear-cut than might at first appear. 

The five studies^''^^'^^'^^'^'^ that failed to find a significant overall rela- 
tion between low WSS and disease are naturally of particular interest, 
since they contradict the current consensus. They all computed flow 
in the geometry that was used to assess disease, and they all used a 
wall thickness-related metric as an index of disease severity. It is tempt- 
ing to discount them on these grounds, but many studies using one or 
other of these methods did find a correlation, as did two studies^^'^^ 
that used both methods. It could equally be claimed that four of the 
five are methodologically the best studies because they used a quantita- 
tive, point-by-point comparison between the two variables (and the 
only other study^^ using point-by-point analysis also failed to obtain a 
correlation with this method, even though a correlation was obtained 
after axial averaging of the data). The fact that no point-by-point analysis 
supported the consensus is rendered even more striking by the fact that 
such methods tend to over-estimate correlations (see below). 

Since this review did not distinguish between different vascular 
beds, it is important that the five studies were not restricted to a 
single-vessel type but included coronary and carotid arteries. 
Equally, since this review did not distinguish between human disease 
and lesions occurring in hypercholesterolaemic animals, with or 
without genetic modification, it is important that all five negative 
studies used human vessels; the investigations which used animal 
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models conformed to the consensus, and that is also true of a study 
not meeting our inclusion criteria, where a flow-modifying cuff was 
placed around the carotid artery of hypercholesterolaemic mice — 
lesions developed in regions of low TAWSS.^°'^'' (Interestingly, 
lesion type appeared to depend on the degree of oscillation of the 
shear.) 

A number of issues will need to be addressed by future studies 
in order to improve and refine the existing evidence. 

5.1 Are we capturing the relevant flow 
metrics? 

Do the conventional metrics — WSS, TAWSS, OSI, and RRT — capture 
physiologically relevant features of the flow? And do they quantify the 
essence of the expression 'disturbed flow' that is so commonly used 
in the biological literature? Or should we explore other parameters, 
such as the directional OSI suggested by Chakraborty et o(.?*'^ 

A number of authors have suggested an indirect link between flow 
and disease that depends on flow affecting fluid-phase mass transport 
of solutes (oxygen, LDL) near the wall, rather than direct effects of 
flow on the wall itself WSS and transport are related, but 
idealized models show that their patterns need not be identical."'^^ 

5.2 What level of data reduction is required 
and justifiable? 

The use of different data reduction techniques to decrease the 
amount of data-processing (e.g. by limiting the analysis to a subset 
of the available data, averaging to combine data, or thresholding) 
complicates the evaluation of the evidence. Authors should make it 
possible to evaluate results after minimal processing, for example, 
by providing figures that allow a visual comparison, even when statis- 
tics are applied only to reduced measures. 

Point-by-point comparisons represent an ideal. Nevertheless, there 
remain limits to the spatial resolution that can be attained. In particu- 
lar, cross-correlation between adjacent data points must be minimized 
to prevent over estimation of statistical significance.^^ Errors in image 
registration are also limiting. 

5.3 Are subject-specific studies introducing 
errors? 

Studies of the effects of WSS on the initiation of atherosclerosis should 
compute flows using disease-free geometries. On the other hand, the 
large effect that small changes in geometry can have on WSS patterns 
suggests that subject-specific investigations are best. Currently, these 
two requirements are hard to reconcile: they require either that 
effects of disease on geometry are negated naturally (by compensa- 
tory remodelling) or removed artificially (by hard-to-validate 'virtual' 
techniques), or that longitudinal studies are carried out over 
periods — at least for human disease — of many years. 

These problems could be resolved by the development of techni- 
ques that allow disease to be mapped at a stage before intimal thick- 
ening occurs, and by methods that could ethically be used in such 
minimally affected subjects. Until such techniques are available, it 
may be more appropriate to drop the subject-specific approach and 
instead use different study groups: a diseased one for the assessment 
of lesion patterns and a disease-free one for the assessment of flow. 
The essential requirement of sufficiently large samples in the haemo- 
dynamic group has been fulfilled in only one study of those reviewed, 
and in none of those that considered human vessels. 



5.4 Have all disease distributions been 
considered? 

Disease distributions sometimes seem inconsistent. This is illustrated 
by the original debate over the high and low shear theories: Fry^ 
observed lesions at the downstream margin of aortic side branches, 
a site that Caro et o(.^ found to be spared of disease. As both theories 
assumed that this region experiences high WSS, the debate con- 
cerned the distribution of lesions rather than characteristics of the 
blood flow. 

Since the high shear theory was predicated on the distribution seen 
in animal models while the low shear theory was based on a study of 
human post-mortem material, the discrepancy was attributed to dif- 
ferences in plasma cholesterol levels.^'^^ However, it has subsequently 
emerged that lesion patterns change with age in rabbit and human 
aortas, and in a parallel way;^' the discrepancy might be explained 
better by such changes. Lesions can also show different patterns in 
one species at a single age. In the mouse, they occur upstream of 
branches originating in the aortic arch, all around branches of the 
thoracic aorta, and downstream of branches in the abdominal 
aorta.^°'^^ 

The different patterns provide an additional challenge for haemo- 
dynamic theories of atherogenesis. The fact that the large majority 
of studies concerning shear and atherosclerosis have concentrated 
on relatively few locations in mature arteries may have led to an 
underestimate of the variety of mechanical factors that can trigger 
disease. 

Supplementary material 

Supplementary material is available at Cardiovascular Research online. 
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